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We study the size and the hfetime distributions of scale-free random branching tree in which k branches are 
generated from a node at each time step with probabihty ^j, ~ k^"^. In particular, we focus on finite-size trees in 
a supercritical phase, where the mean branching number C = Y,kkqk is larger than 1. The tree-size distribution 
p{s) exhibits a crossover behavior when 2 < y < 3; A characteristic tree size Sc exists such that for s <C Sc, 
p{s) ~ j^T/(T^') and for i 2> ic, p{s) ~ s^'^l'^ exp(— i-/.Sc), where Sc scales as ~ (C— 1)^(1'^ ')/(T^2) pgj- y > 3^ 
it follows the conventional mean-field solution, p{s) ~ i^-'/^ exp(— i/iv) with Sc ^ {C — 1)^^. The lifetime 
distribution is also derived. It behaves as i{t) ~ f^(T^')/(l'^2) for 2 < Y< 3, and ~ f^^ for Y> 3 when branching 
step r<Cfc~(C— 1)^', and l{t) ~ exp(— ?/fc) for all Y> 2 when t^tc- The analytic solutions are corroborated 
by numerical results. 



I. INTRODUCTION 

A tree is a graph with no loop within it. Owing to the sim- 
plicity of its structure and amenability of analytic studies, tree 
graph has drawn considerable attentions in many disciplines 
of scientific researches. Scale-free (SF) random branching 
tree, in which the number of branches k generated from a node 
is stochastic following a power-law distribution, qi^ k^"^, is 
particularly interesting here. Such trees can be found in var- 
ious phenomena such as the trajectories of cascading failure 
in the sandpile model on SF networks fT], epidemic spread- 
ing on SF networks [H, aftershock propagation in earth- 
quake lid, random spanning tree or skeleton of SF net- 
works f&\, phylogenetic tree fT], etc. Here, SF network is 
the network with the degree distribution following a power 
law P,i{k) ~ /k"'' S H [mi. So far, several analytic studies 
have been performed to understand structural properties of SF 
branching trees 1 11]. However, most works are focused on the 
critical case, where the mean branching number C = Y,k is 
equal to 1, motivated by universal feature of scale invariance 
observed in nature and society. 

Recent studies, however, show that the structure of real- 
world networks may have been designed upon supercritical 
trees [6]. Supercritical trees, where the mean branching num- 
ber C > 1, turn out to act as a skeleton of some fractal net- 
works such as the world-wide web. Here skeleton [12] is 
defined as a spanning tree forme d by edges with highest be- 
tweenness centrality or loads |T3','l4']. A supercritical branch- 
ing tree can grow indefinitely with a nonzero probability, 
which is the most marked difference from critical (C — I) or 
subcritical (C < 1) tree that cannot grow infinitely. Moreover, 
the total number of offsprings 5(f) generated from a single root 
(ancestor) up to a given generation t can increase exponen- 
tially in supercritical trees and this is reminiscent of the small- 
world behavior: The mean distance between nodes scales log- 
arithmically as a function of the total number of nodes ill 1(1 . 

Due to the mean branching number being larger than 1, 
some supercritical trees may be alive in a very long time limit. 
The tree-size distribution of those surviving trees in the su- 
percritical phase has been derived in the mean-field frame- 



work lITsll . which follows a power law, p{s) s^^. Here, we 
consider finite-size trees in the supercritical phase. In spite 
of the large mean branching number, some trees do not grow 
infinitely even in the supercritical phase. For such finite-size 
trees in the supercritical phase, we derive the tree-size and 
the lifetime distributions using the generating function tech- 
nique l(l6ll . Distinguished from the critical case, the generat- 
ing function of the tree-size distribution exhibits two singular 
behaviors in the supercritical phase and thereby a crossover 
behavior of the tree-size distribution can arise when 2 < y < 3. 
We present in detail the derivation of all these analytic so- 
lutions in the following sections. The tree-size and lifetime 
distributions predicted by analytic solutions are confirmed by 
numerical simulations. This is important in itself for under- 
standing the branching trees whose structure changes drasti- 
cally depending on the phase. Since the branching tree ap- 
proach can be applied to numerous systems, our results should 
be useful for future diverse applications as well. 



II. TREE-SIZE DISTRIBUTION 

Let us consider the branching process that each node gen- 
erates k offsprings with probability qk. 



1 



qk 



for k = Q, 



■k-^ for k>\, 



(1) 



where C is constant in the range of < C < ^(y — 1)/^(y) with 
the Riemann-zeta function ^(x), and y is larger than 2, ensur- 
ing that ^(y — 1) is finite. Then, C is automatically identical to 
the mean branching number, i.e. the average number of off- 
springs C ~ YX=o^qk generated from a node. When C < 1, 
the number of offsprings decreases on average as branching 
proceeds and it vanishes eventually. Thus, branching tree has 
finite lifetime with probability one. When C > 1, as branching 
proceeds, the number of offsprings can increase exponentially 
with non-zero probability. The case of C = 1 is marginal: Off- 
springs persist, neither disappear nor flourish on average. A 
branching tree generated through the stochastic process (1) is 
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a SF branching tree, because its degree distribution follows a 
power law, Pci{kd) ~ asymptotically. Degree of each 
node in the tree is related to the branching number k of that 
node as A:^/ = + 1 but for the root, kd = k. 



A. Generating function method 

A tree grows as each of the youngest nodes generates their 
offsprings following the probability (7^. in Eq. ([T]i- This evo- 
lution is regarded as a process in a unit time step. When a 
node generates no offspring with probability qo, it remains in- 
active in further time steps. We define pt{s) as the fraction 
of trees with total number of nodes s at time t. By definition, 
p{.){s) = 5i,i. Then, pt+\ [s) can be written in terms of pt{s) as 

k=Q si,S2 it 

(2) 

Defining the generating functions, Q,(co) = T^^olk^^ and 
'Pt{y) — T.7=i Pt{s)y\ and applying them to dU, one can ob- 
tain that 

i'r+i(y)^y(l{My))- (3) 

Let us consider the tree-size distribution in the t ^ °o limit, 
i.e., p{s) = lim,^oo pt{s) and its generating function (P{y) = 
lim(^„o ^t{y)- However, some trees may grow infinitely in the 
supercritical phase, which makes £P (y) = ill-defined 
at y = 1 . So we limit the summation in fP (y) over finite trees 
only, i.e., = Y^hnite sP{^)y'''- This is equivalent to defining 
2'(1) — limy^i£P(y). Then, Eq.Q gives the relation in the 
f ^ 00 limit, 

¥{y)^y(l{¥{y)). (4) 

The next step is to extract a singular part of rP (y) from Eq. (|4|l, 
and then to derive the behavior of p{s) for 5^1. 

The power-law form of in Eq. ([U results in the expansion 
of Q, (co) around co 1 : 



Q.((0) = 1-C(1-C0) 

A(y)(1-co)^-' 



B{y) 



(l-co)2 + ... 
(Y / integer) 



(5) 



^(l-co)T iln(l-co) (y = integer) 

where B{y) = C[C,{y - 2)/i;{y - 1) - 1], and A{y) = Cr(l - 
Y)/C(Y~ 1 ) with the Gamma function r(x). The inverse func- 
tion y — ¥^^{(0) is then expanded as 



CO 



Q,(«) 
A(y)(1-co)T^-i 



l+A(l-(o)-^(l-(o)2^ 
(Y ^ integer) 



i^(l-co)T-iln(l-co) (Y = integer) 



(6) 



where A = C — 1 . We recall that A is positive (negative) in 
the supercritical (subcritical) regime and in the critical case. 
Here we focus on the supercritical case of A > and being 
very small, but the obtained result can be naturally extended 
to large-A cases. 








CO* 



CO 



FIG. 1: (Color online) Schematic plot of the function y = co/q,(co) 
in the supercritical phase. The dy/da) = occurs at O) = CO* < 1. 



B. The singularity at y = y* > 1 

Let us investigate how y behaves as CO decreases from 1 to 
0. For A > 0, as CO decreases from 1 to CO*, y increases from 
1 to y* and then decreases to zero as shown in FiglT] where 
CO* satisfying (ii/c/co)[co/Q.(co)]|(o=(o* = locates less than I. 
This feature is distinguished from the solution CO* = 1 for the 
critical case. It is obtained that CO* depends on A as 

{A for Y>3, 

A/ln(l/A) for Y=3, (7) 
Ai/(Y-2) for 2<Y<3. 

The value y*, determined by the relation y* = co*/q,(co*), lo- 
cates at 

A2 



for Y>3, 



y*-l=5*-<^ AVln(l/A) for Y=3, 



(8) 



^(Y-l)/(Y-2) foj. 2<Y<3. 



The curve y = co/q,(co) in the region CO > CO* is just the ana- 
lytic continuation of the inverse function y = fP^'(co) that is 
analytic for CO < CO* iITtIi . 

The right-hand-side of Eq. (|6]l for co < CO* is expanded 
around CO* as 



when CO is close to co* such that 

A,+i(y) 



max ■ 



(co* - co) < L 



(9) 



(10) 



«>2 D„(Y)(n+l) 
Here D„(y) is the nth derivative of co/q,(co) at CO*. For « = 2, 
-B(y) for Y>3, 

DiXi) ^ { InA for Y=3, (11) 

_A(T^-3)/(y-2) for 2<Y<3. 
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This result is used for future discussions. Keeping only the 
quadratic term (co^, — co)^ in Eq.®, one obtains the leading 
singular behavior of rP {y) at y,, 



CO = !P (y) 



' CO* — - 



(12) 



In fact such a square-root singularity at y — y.^ is generic 
regardless of the form of the branching probability when 
10 + 11 < 1 ITtIi . yielding the asymptotic behavior of p{s) 
given by 

pis)^b{A)s-^^^exp{-s/s,), (13) 

where the coefficient b{A) - a-(t-3)/P(T-2)1 for 2 < y < 3, 

1 / ^ln(l/A) for y = 3 and constant for y > 3, and 5* — 
{lny,)-\ 

C. The singularity at y = 1 

When CO is far from CO* such that the linear term with the 
coefficient A is not comparable to the next-order term, an- 
other singularity becomes dominant. The next-order term is 
the quadratic term for y > 3 and the non-analytic term for 

2 < y < 3. To be precise, if the condition, 1 — CO ^ A for y > 3, 
-(1 -co)ln(l - co) > Afory=3, and 1 - co > A'/^l'-^) for 
2 < y < 3, holds, then the linear term is negligible compared 
with the next order terms, and then Eq.® is reduced to 



^(l-co)2 for y>3, 

-co)2ln(l -co) for y=3, 
A(y)(l -co)T'-i for 2<y<3. 

The generating function £P (y) then behaves as 



(14) 



CO = rp (y) ~ 1 



2(i-y) 

B(Y) 

4(1 -y) 

|ln(l-v)l 

^\ 1/(7-1) 
My) 



for y>3, 
for y=3, 
for 2<y<3. 



(15) 



From this result, one can obtain the tree-size distribution as 

s^^/^ for y>3, 

p{s)^{ i-3/2(lni)-i/2 for y^3, (16) 
,-T/(Y-i) for 2<y<3. 



D. Crossover behavior between the two singularities 

The two singular behaviors of fp (y) in the forms of Eqs. ( fT2] i 
and ( fTSI ) occurring at y = y* and y = 1 , respectively, enables 
us to determine the ranges of size s where the formulae of 
Eqs.(fT3]) and ( fTSI l are valid. In particular, when 2 < y < 3, the 
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FIG. 2: (Color online) The tree-size distribution p{s) for y = 3.3 for 
various values of C in the scaling form Eq. J17t . Dashed line is guide- 
line with a slope of —3/2. Inset: Dependence of the characteristic 
size Sc on the mean branching number C. 



asymptotic behaviors in Eqs.dTsTi and (fTST i differs from each 
other and thus there should be a crossover behavior in the tree- 
size distribution. 

The ranges of CO in which Eqs. (|9]l and ( fT4b are valid are 
closely related to those of y for Eqs. (fT2l l and (fTSl l and that 
of s for Eqs. ( fT3] l and ( fT6l l. respectively. Here we find those 
ranges of CO, y, and s, and then determine the crossover in the 
tree-size distribution p{s). 

First, we study valid ranges of Eqs. (|9]l, (fT2] |. and (l3[ . 
The coefficient I>„(y) in Eq. ^ behaves as ^ (1 — co*)''^'"" 
for n > y — 1 due to the non-analytic term (1 — co)^^' in 
Eq. ^ when y is not integer. Then, it follows that [Z)„ (y) (n + 
l)/Z)„+i(y)] 1/(1 -co*) = 1/e*. Thus, the condition ^ 
can be rewritten as CO* — <g; CO < co*, where £^ ^ A for 
y > 3, e> - A/ln(I/A) for y = 3, and e> - A^iy-^) for 
2 < y < 3 from Eq.©. The corresponding range of y is 
J* — 5^ <C y < y*, where 5^ is given by ^ A^ for y > 3, 
~ A2/ln(l /A) for y = 3, and - A^^-i'/t^-^) for 2 < y < 3 by 
using Eqs.® and (fT2] i. 

To find valid range of s for p{s) in Eq. (fT3l) . we use the fact 
that the singular functional behavior of CP (y) around y = y is 
determined by that of p{s) around s = S, where y and s are re- 
lated as y" r-^ 1 . Then, one can find that = | ln(3'* — 5;? ) | ^ ^ ^ 
(5* - 5> )" ' , so that sf - A-2 for y > 3, A^^ in( 1 /A) for y = 3, 
and A"(1'"')/(y-2) for2<y<3. For the range the for- 

mula ( fTST l is valid. 

Second, we check the validities of Eqs. (fT4l) . ( fTSI ), and 
( fT6b . Comparing the magnitude of the linear term and the 
next-order term in Eq. we find that Eq. ( fT4b is valid for 
CO ^ 1 — e^*?, where behaves as A for y > 3, A/ln(l/A) for 
y = 3, and A^/'''^^' for 2 < y < 3. The corresponding range 
of y for Eq. ( fTSb is given as y <C 1—5^, where 5^ ~ A^ for 
y > 3, 5< - AVln(l/A) for y = 3, and 5< - a(^-i)/(^-2) ^j. 
2 < y < 3. The corresponding range of s for Eq. ( fTSI ) is s <C sf 
with = |ln(l - 5<)|-i ~ (5<)^' given by s< - A'^ for 



y<3,^<^ 
2<y<3. 



. A-2ln(l/A) for y = 3, and sf - A-(^-1)/(t-2) for 




FIG. 3: (Color online) The tree-size distribution p{s) for y — 2.5 for 
various values of C in the scaling form Eq. ( 119b . Dashed line is guide- 
line with a slope of —5/3. Inset: Dependence of the characteristic 
size Sc on the mean branching number C. 




As already noticed, the crossover sizes sf, sf, and 
are consistent for all values of y within A-dependence, and 
thereby, we use the notation for all of them. The overall be- 
havior of the tree-size distribution is obtained by combining 
Eqs. (T3[ and ( fT6b . For y > 3, there is no need to introduce a 
crossover Thus, it leads to 



7is)r^s-^/^exp{-s/sc) (y>3), 



(17) 



for all s. And Sr 



As A increases, the cut-off Sc de- 



creases and the exponential-decaying pattern prevails. 
When Y= 3, p{s) is given by 



i"^/2(lni)"'/2 for s^Sc, 



,-3/2 



exp(-i/if) for s:$>Sc, 



(y=3) (18) 



where Sc ^ A ^ ln(l/A). Similarly, for 2 < y < 3, we find that 



fors < Sc, 



,-3/2 



exp{~s/sc) for s:^Sc, 



(2<Y<3) (19) 



where both Sc ~ A^'^^''/'"'^^). As A ^ 0, Sc diverges, and the 
power-law behavior prevails. 

We invoke numerical simulations to confirm our analytic 
solutions. Figs. |2] and [3] show the tree size distributions for 
Y= 3.3 and Y= 2.5 in the scaling forms, Eq. il% and Eq. ( fT9l ). 
respectively. The data are well collapsed into the predicted 
formulas for different C values for both cases. 



III. LIFETIME DISTRIBUTION 

Next we solve the lifetime distribution £{t). This is defined 
as the probability that the branching process stops at f . To de- 
rive £{t), we first introduce the probability that the branching 
process stops at or prior to time f, denoted by r{t). Then £{t) 



FIG. 4: (Color online) Schematic plot of the function dr/dt = 
Q,(r) — r in the supercritical phase. {d/dr){dr/dt) = occurs at 
at which dr/dt is denoted as rj.. 



is given as £{t) = r(f + 1) — r(t). The probability distribution 
r{t) is related to r{t — 1) as 

rit)^Y.^kHt-l)]'=Clir{t-l)). (20) 

*:=0 

Thus, we are given approximately a differential equation for 



dr{t) 
dt 



:^(f) = (l(r(f))-r(f), 



(21) 



Expanding the right hand side of Eq. (|2TI) around r = 1, one 
can see its asymptotic behavior. Using Eq. (|5]l again, we find 
dr/dt in the long time limit as follows: 

dr , , , , B(y) , 

- - (i(,)-,= -A(l-r) + ^(l-r)2 + ... 



A(Y)(l-r)^-l 

(-i)T 



(Y 7^ integer) 



r7y-(l-r)^-iln(l-r) (y- integer) 



•(22) 



What we can see in this relation is that the value of r' is zero 
at r = 1 . It decreases as r decreases until it reaches where 
{d/dr)\Ci(r) — r\\r=r, — holds. Passing r*, r' increases as r 
decreases further, crossing the r' = as shown in Fig.lD 

First, as in the case of a)/Q,(co), two singularities exist in 
qSj) — f- For r close to r,, Eq. (l22l) is expanded as 



n=2 



(23) 



where = Q,(r*) — < and G„(y) is the n-th derivative of 
Q,(r) — r at r* . When r is close to such that 

G„+i(y) 



max„>2 



g„(y)(« + i) 



(r*-r)<l, 



(24) 
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one may neglect higher order terms, keeping only the 
quadratic term in — r as 



The solution to the above differential equation is 



2a 



r(t) ~ r(oo) ■ 



(25) 



(26) 



where r{oo) = r^, — a and a — \Jiy^\lG2 (y), and — 
l/^2|r;|G2(Y)- The Ufetime distribution ^(f) = r'(f ) is then 
given by 



(27) 



Second, following the same steps taken for the singularities 
of T {y), we find another approximate relation between r' and 
r in the region of r{t) where the next order term in Eq. ( |22] | is 
much larger than its linear term as follows: 



( ^(1-^)2 for Y>3, 

^~<^ -i(l-r)2ln(l-r) for y=3, 

I A(y)(l -r)T-i for 2 <y< 3. 

Their solutions are, in long time hmit, given by 

f"' for y>3, 

\~r{t)^ { f-'(lnf)-i for y=3, 

f-i/(T-2) for 2<y<3. 



(28) 



(29) 



From these results, the lifetime distributions are obtained as 

{r^ for y>3, 

f-2(inf)-i for y=3, (30) 
j-(Y-i)/(Y-2) foj. 2<y<3. 

Different behaviors of the lifetime distribution shown in 
Eqs. dZTl l and (|30] | suggest the presence of a crossover behav- 
ior The characteristic time that distinguishes the two behav- 
iors for given y can be found by considering the valid ranges 
of t for Eqs. (|27] | and (l30l l. respectively. When the condi- 
tion of Eq. ( l24l i is fulfilled, Eqs. ( |26] | and dZTl i are valid. The 
condition is approximately represented in different form of 

- r < 1 - since G„(y) - (1 - r*)^-'-". From Eq. (|22]l, 
one can find the value of 1 — for different y's: 1 — ~ A 
for y > 3, A/ ln(l /A) for y = 3, and Ai/(l'-2) for 2 < y < 3, re- 
spectively. Applying these conditions to Eq. (l26T l. it is found 
that Eqs. (|26] | and ( l27l) are valid if f :§> f*i with f*i ~ A^' irre- 
spective of y as long as y > 2. 

Eqs. ( |29] l and (|30] | are valid when the linear term is much 
smaller than the next order term, which is satisfied when 1 — 
r > A for y > 3, 1 - r > A/ln(l/A) for y = 3, and 1 - r > 
l^^liy-^) for 2 < y < 3, respectively. Applying these conditions 
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FIG. 5: (Color online) The lifetime distribution i'(f) for Y= 3.3 in the 
scaling form Eqs. l l31t and l l32t . Dashed line is guideline with slope 
—2. Data for small t are deviated from the data collapse, indicating 
that our solution is valid for large t only. Inset: Dependence of the 
characteristic time tc on the mean branching number C. 
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FIG. 6: (Color online) The lifetime distribution l{t) for y = 2.5 in the 
scaling form Eqs. l |31| l and l |32| l. Dashed line is guideline with slope 
—3. Data for small t are deviated from the data collapse, indicating 
that our solution is valid for large t only. Inset: Dependence of the 
characteristic time tc on the mean branching number C. 



to Eq. (|29] l leads commonly Xo t <^tt2 ^ ■ One can find 
that the two characteristic times f*i and f*2, and scale in 
the same manner, so that they are denoted as tc commonly. 
Therefore, we conclude that the lifetime distribution behaves 
as 

{t^^ for y>3, 

f-2(inf)-i for y=3, (31) 
j-(Y-i)/(T-2) for 2<y<3, 

when t <^tc ^ Ar^ , and 

l{t)^e-'l'' for y>2 (32) 

when t tc- The analytic solutions for the lifetime distribu- 
tion are checked by numerical simulations in Figs. |5]and|6] 
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Data in small t regime are somewhat deviated from the data- 
collapsed formula, indicating that our solution is valid in large 
t regime. 

IV. CONCLUSIONS AND DISCUSSION 

Our main results are Eqs. dTTl l. dTsl l. and ( fT9] l for the 
tree-size distribution when trees are finite: Contrary to the 
case of Y > 3 for which the tree-size distribution p{s) be- 
haves as ~ ,s'^-'/^exp(— i/if) for all s with Sc ^ {C — 1)^^, 
a crossover behavior occurs at Sc ~ (C — 1)^(t-i)/(Y-2) ^^j. 
2 < y < 3. For 5 < Sc, p{s) ~ s-'l^/(T-i) and for s > Sc, 
p{s) ^ s^^l^ e.x^{—s / Sc). This result is complementary to 
the previous mean-field solution pinf (s) ~ s^^ for infinite-size 
tree. From our solutions, it is noteworthy that the characteris- 



tic size Sc increases as the exponent y approaches 2. This leads 
to an interesting result: A larger-size tree can be generated for 
smaller value of the exponent y. However, the probability to 
have such a large-size tree becomes smaller as the exponent y 
approaches 2, because the exponent y/ (y — 1 ) for the tree-size 
distribution p{s) becomes larger. 

The lifetime distribution also exhibits a crossover behavior 
at f,. - (C- 1)"*. It follows Eq. (EB for f < r,. and (O for 

t > tc. 
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